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Using selection rules imposed by the Pauli principle, we classify pairing 
correlations according to their symmetry properties with respect to spin, mo- 
mentum, and energy. We observe that inhomogeneity always leads to mixing 
of even- and odd-energy pairing components. We investigate the supercon- 
ducting pairing correlations present near interfaces between superconductors 
and ferromagnets, with focus on clean systems consisting of singlet supercon- 
ductors and either weak or half-metallic ferromagnets. Spin-active scattering 
in the interface region induces all of the possible symmetry components. In 
particular, the long-range equal-spin pairing correlations have odd-frequency 
s-wave and even-frequency p-wave components of comparable magnitudes. 
We also analyze the Josephson current through a half-metal. We find an- 
alytic expressions and a universality in the temperature dependence of the 
critical current in the tunneling limit. 

PACS numbers: 7445,+c, 74.20.Rp, 74.50. +r 

1. INTRODUCTION 

The rich physics specific for boundary regions between superconduct- 
ing and ferromagnetic materials has recently been probed in a series of 
e^penments^MMMMimmmmmmmmmimmm Apart from confirm- 
ing earlier theoretical predictions^ 2 ^^^ they provide deep insight into 
the coexistence of the two types of order and have inspired new ideas in the 
emerging field of spin electronics. In the boundary region the characteristic 
correlations known from the proximity effect in normal metals 2 ^ are induced, 
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but in addition to the usual decay, they show an oscillating behavior.^ The 
length scales for decay and oscillations are set by the magnetic length.^ It 
is smaller than the decay length in a normal metal, for ballistic systems by 
the ratio between temperature and exchange energy, ksT/J, and for dif- 
fusive systems by the factor ^ksT/J. This leads to modifications of the 
density of states^^ZEUGS] an( j Josephson effect through ferromagnets 
For usual ferromagnets J is considerably larger than fc^T, and the proximity 
effect is hard to observe. A breakthrough came with the advance of dilute 
ferromagnetic alloys with rather low spin polarization, such as Pdi_a;Ni x 
or Cui-zNiz, where magnetic Ni ions are integrated into a non-magnetic 
matrix. A high level of control has been reached with heterostructures con- 
taining such "weak" ferromagnets (low spin polarization). E.g., it became 
possible to spatially resolve properties on the scale of the magnetic length 
and to observe the proximity effect 

In parallel, it remained an important goal for the further development of 
spin electronics to create and investigate heterostructures^ ^^ ^M^ where 
"strong" ferromagnets with large exchange splitting of the bands and high 
spin polarization are in contact with superconductors. In the extreme case, 
for so-called half-metallic materials, the polarization reaches values close to 
100 %. Half metals are metallic in one spin direction with respect to a cer- 
tain spin quantization axis, and semiconducting or insulating in the other. 
In the course of this research a so-called "long-range proximity effect" was 
discovered JESIHEl CTI^ E^IMlEll w hich is governed by a length scale typical of 
the proximity effect in a normal metal, rather than the magnetic length of 
the ferromagnet involved. Such long-range proximity amplitudes have their 
origin in the mixing between singlet and triplet pair amplitudes inherent 
to pairing in the presence of broken spin-rotation symmetryP^ElGIS] Triplet 
correlations are classified according to their projection on the spin quantiza- 
tion axis that renders the quasiparticle bands in the ferromagnet diagonal. 
There are three corresponding amplitudes, m = 0, dbl. Of those the equal 
spin pair amplitudes, m = ±1, lead to the long-range proximity effect, as 
pairing occurs within the same spin band and is unaffected by the large ex- 
change energy J. This explains the penetration of triplet pair amplitudes 
into the ferromagnet, once they are created G23E2E21 However, the questions 
how they are created in the first place and what are their magnitudes remain 
active research topics \MEMMEE!MMiMMMMI] 

The creation mechanisms for pairing correlations inside the ferromagnet 
differ in the ballistic limit significantly for weak and strong ferromagnets. For 
the following discussion we will concentrate on singlet superconductors. In 
the case of a weak ferromagnet the singlet superconductor induces in the first 
place a mixture of singlet and m = triplet pairs in the ferromagnet. Both 
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penetrate only on the short magnetic length scale, but various mechanisms 
may create long-range m = ±1 triplet components. Examples are (i) a 
magnetic domain wall near the interface within a distance of the order of the 
magnetic length EES or ^ ^ w0 ferromagnets with misaligned quantization 
axes separated by a singlet superconductor with a thickness of the order 
of the superconducting coherence length^SKD or ( m ) a spin-active interface 
that allows for spin-flip processes.^ In all cases, magnetic inhomogeneities 
mix the triplet pair components and create the long-range equal-spin pair 
amplitudes. 

In the case of a strong ferromagnet the roles of the ferromagnet and the 
superconductor are reversed. Here, in the first place the ferromagnet acts as 
a source for spin-polarization of Cooper pairs in the superconductor. This re- 
sults in a boundary layer with coexisting singlet and triplet amplitudes near 
the interface extending about a coherence length into the superconductor. 
The important mechanisms here are the spin-mixing termiP^ESHIMIMSK^ 
(often called spin-rotation) in the reflection and transmission amplitudes of 
the surface scattering matrix. Such spin-mixing effects arise as a conse- 
quence of the different matching conditions for spin-up and spin-down wave 
functions at the interface.^ Consequently, the creation of triplet pair am- 
plitudes is entirely the result of the interface properties, taking place in an 
interface region that is of similar size as the magnetic length. Spin-mixing is 
most effective at interfaces with strong ferromagnets, increasing in strength 
with growing spin-polarization of the ferromagnet. Long-range triplet com- 
ponents are created when spin-flip centers are present in the interface region. 
This mechanism even works in the presence of completely polarized ferro- 
magnets, since the triplet correlations are created entirely within the super- 
conductor, and only after their creation penetrate into the ferromagnet 
In addition, the magnitude of the triplet correlations at the interface is pro- 
portional to that of the singlet amplitude at the interface, and both are 
insensitive to impurity scattering (in contrast to the decay behavior away 
from the interface).^ 

For intermediate spin polarizations the two creation mechanisms for 
triplet pairing, with strengths depending on microscopic details of the in- 
terface and domain wall structures in the ferromagnet, compete and are 
difficult to characterize theoretically. However, the two limiting cases of 
weak and strong spin polarization can be treated within a controlled ap- 
proximation, namely the quasiclassical theory of superconductivity.^"^ It 
relies on a separation of two energy scales, the superconducting gap A and 
the Fermi energy Ef, and can be extended to superconductor-ferromagnet 
heterostructures when the exchange energy J lies either in the low- or in 
the high-energy range, corresponding to weak and strong ferromagnets, re- 
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spectively. For weak ferromagnets, this means that Fermi surface properties 
like the density of states at the Fermi level in the normal state, Nj, and the 
Fermi velocity, Vf, remain unchanged to lowest order in the small quasiclas- 
sical expansion parameter J/Et; note that the magnetization in the normal 
state is a first order term in this parameter and is given by 2/j.BNfJ, where 
HB is the Bohr magneton. On the other hand, for strong ferromagnets all 
effective interactions and indeed the quasiparticle band structure itself will 
be strongly modified by the presence of the exchange splitting. It is not 
possible a priory to describe the crossover between the two limits within 
quasiclassical theory. 

The outline of this article is as follows. In Sec. 2 we describe how the 
Pauli principle leads naturally to a classification of superconducting corre- 
lations according to their symmetries with respect to spin, momentum, and 
energy. We then show by two examples that all these types of correlations 
are indeed induced at superconductor- ferromagnet interfaces. In Sec. 3.1 
we analyze the case of a weak ferromagnet in contact with a singlet su- 
perconductor through a spin-active interface barrier. We also illustrate the 
difference between the spatial dependences of short-range and long-range 
proximity amplitudes. In Sec. 3.2 we consider the case of a strong ferromag- 
net. In particular, we show results for the proximity effect between a singlet 
superconductor and a half-metal, and for a Josephson junction involving a 
half metal. We show that an analytic treatment is possible for the case of 
small transmissions and a small spin-mixing angle. We find that in these 
limits, a previously discovered^ low-temperature anomaly in the critical 
Josephson current is independent of the interface parameters, and is a direct 
consequence of the different symmetry properties of the pairing correlations 
compared to the case of a normal metal between two superconductors. 

2. SYMMETRY CLASSIFICATIONS 

There has been considerable work that formulated the physics of pair- 
ing correlations in diffusive ferromagnets in terms of one particular pairing 
component, that has odd- frequency, s-wave, equal-spin triplet symmetry.^ 
The question arises to which extend this component is important for sys- 
tems with weak or intermediate impurity scattering. We will show in the 
following chapters that in fact four different symmetry components exist in 
ballistic heterostructures, and that they are comparable in size. We start 
with a general classification of pairing correlations. 

Superconducting correlations are quantified by the anomalous Green's 
function 

^9 (ri , 7i ;r 2 ,T 2 ) = (T r * a (ri, ri)^(r 2 , r 2 )> . (1) 
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It is a matrix in spin space and depends on two coordinates and, in the 
Matsubara technique, on two imaginary times. The Pauli principle requires 
that this function changes sign when the two particles are interchanged, 

J r ai3(ri,T 1 ;r 2 ,T 2 ) = -J r (3a (f 2 ,T 2 ;f 1 ,T 1 ). (2) 

This well-known condition follows directly from Eq. (JTJ) and the anti- 
commutation relations for the field operators. For homogeneous systems 
T depends only on relative coordinates r = r\ — r 2 and r = t\ — t 2 , i.e. it 
follows F a f){r,T) = — Fpai— r, — t), and after Fourier transformations 

Fa/3(P, £n) = -T{3 a {-V, ~^n)- (3) 

For inhomogeneous systems this equation holds for each set of center coor- 
dinates. The symmetry restriction Eq. ([3]) in spin, momentum p, and Mat- 
subara frequency, e n = (2n + l)vrT, can be satisfied in four different ways, 
listed in Table [TJ Analytical continuation to the complex z-plane leads to 
TufiiPiZ) = —Fj3 a (—p,—z); in particular, retarded and advanced functions 
are related as T^p, e) = — Fp a (— P, — e)P2 

The spin part of Eq. ([3]) can be divided up into singlet and triplet sectors 

Fapip, e n) = Faff, {^y) a j3 + ^t(p, Cn) " {^^y) a p , ( 4 ) 

where a = (a x , cr y , a z ) is the vector of the three Pauli matrices. The singlet 
spin matrix {ia y ) a p is odd under the interchange while the three 

triplet matrices (ai<Ty) a s are even. Some insight about the symmetries in 
the momentum- and frequency-domains can be gained by considering the 
equal-time correlator EMMM] E.g., for the spin-triplet case the Pauli principle 
imposes that the following sum must vanish 

F t (?= 0, r = 0) = TJ2 E e «) = °" ( 5 ) 

When the orbital part is odd the electrons avoid each other in real space, 
the equal-time correlator F t (p,r = 0) = T^ tn J : t{p-,£n) can be finite, and 
the correlator 0~t{p,€n) is even in frequency. On the other hand, when the 
orbital part is even, the correlator J^tip, £n) is odd in frequency, and electrons 
avoid each other in time. 

We summarize the symmetry classes in Table [TJ The usual spin singlet 
s-wave orbital symmetry in a BCS superconductor is of type A, while the spin 
triplet p-wave orbital symmetry superfluid formed in 3 He is of type C. Type 
D was first considered by BerezinskiP^ in connection with early research on 
superfluid 3 He. Finally, type B was considered in connection with uncon- 
ventional superconductors by Balatsky, Abrahams and others E2E2I54J155I 
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Table 1 

The Pauli principle, requiring the pair correlation function to be an odd function 
under the exchange of two electrons, can be met by properties in spin, orbital, or 
frequency space as follows: 



type 


spin 


momentum 


frequency 


overall symmetry 


A 


singlet (odd) 


even 


even 


odd 


B 


singlet (odd) 


odd 


odd 


odd 


C 


triplet (even) 


odd 


even 


odd 


D 


triplet (even) 


even 


odd 


odd 



So far we have only considered the correlation function. In order to 
obtain the gap function, additional knowledge of the pairing interaction 
\(p,p', Sn, e' n ) is required (here for the singlet case) 

A s (p,e n ) = T^^Xsip,^' ,e n ,^ n )T a {p' : e' n ). (6) 
< V 

Clearly, the symmetry of the pairing interaction dictates the symmetry of 
the gap function through the projection of the correlation function in the 
gap equation. In particular, if the superconducting correlations are odd in 
frequency a frequency-dependent pairing interaction (due to strong retarda- 
tion effects) is needed to obtain a non- vanishing gap and a superconducting 
transition. 

3. TRIPLET PAIRING IN CLEAN S/F STRUCTURES 

3.1. Weak Ferromagnets 

To illustrate how triplet correlations are induced we study a 
superconductor-weak ferromagnet heterostructure in the ballistic transport 
regime and include spin-active interface scattering. For simplicity we con- 
sider temperatures near the superconducting critical temperature; however, 
the main results of this section require only the pairing amplitudes in the 
ferromagnet to be small, and apply with minor modifications to any tem- 
perature. In this case, within quasiclassical approximation, the anomalous 
Green's function follows from the (linearized) Eilenberger equations S51HS1 

(v r V + 2e n )f s = 2vrAsgn(e n )-2iJ/f 2 , (7) 
{v f ■ V + 2e n ) f u = -2iJf s , (8) 
(v r V + 2e n )f t± = 0. (9) 

Here, Vf denotes the Fermi velocity for the quasiparticles in the respective 
material. The superconducting gap A is non-zero only in the superconduc- 
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Fig. 1. (a) The spin quantization axis jl of the interface spin-active scattering 
matrix forms an angle a with the exchange field in the ferromagnet and has a 
perpendicular component ji± . (b) The direction of the momentum pf is expressed 
in spherical coordinates through the angles 9 P and <p p . 



tor, while the exchange field J is non-zero only in the ferromagnet. We 
assume that the ferromagnet has a single domain and use the direction of 
its exchange field Jz as spin quantization axis. The spin- active interface can 
have a different spin-quantization axis, that we denote ft (see below). 

The main features of the set of differential equations ([7])- ([9]) are: (i) 
The inhomogeneity of the equation for ft z requires both J and f s to be 
present. Thus, this component naturally emerges in a ferromagnet coupled 
to a singlet superconductor, (ii) The eigenvalues of the f s -ftz sub-system 
for a particular Vf are given by h~ = 2(|e n | ± iJ)/vf. Thus, both the 
singlet f s and the triplet ft z oscillate on the clean-limit magnetic length 
scale £j = Vf/2J, and decay exponentially on the length scale £ n = Vf/2\e n \; 
the latter is dominated by the lowest Matsubara frequency, eo = vrT, and 
occurs on the clean-limit normal-metal coherence length scale £t = Vf/2irT. 
(iii) The equation for f t ± = (ftxifty) is decoupled from the others and is 
homogeneous. Therefore, the presence of these components requires spin- 
active interface scattering, (iv) The equations for ft± do not contain the 
exchange field, and these components are monotonic decaying functions on 
the scale £„. 

The Eilenberger equations are solved by integrating along trajectories 
Vf(j)f) with an initial condition at the starting point of the trajectory. We 
consider the three-dimensional case and introduce a spherical coordinate 
system as shown in Fig.[T](b). The superconductor occupies the region x < 0, 
and the ferromagnet the region x > 0. We assume rotational invariance 
around x, i.e. all quantities will be independent of the variable (ft p . 

For positive Matsubara frequencies, the initial condition for the trajecto- 
ries 6p £ [0,7r/2] deep in the superconductor is f(x — > — oo) = (TrA/\e n \)ia y . 
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For trajectories 9 P S [tt/2,7t] we start the integration deep in the ferromagnet 
with initial condition f(x — > +00) = 0. For negative Matsubara frequen- 
cies, the stable direction of integration is the opposite. After applying the 
boundary conditions described below, the integration proceeds away from 
the interface. 

The boundary conditions at the superconductor-ferromagnet interface 
for the 2x2 spin matrix Green's functions coincide for positive Matsubara 
frequencies with those for the retarded functions.^ Near T c , they read 

fl,mit = <S22/2,in£>22 + S2lfl,in$-21i (10) 

where the indices 1 and 2 refer to the superconducting and ferromagnetic 
sides, respectively, and 'in' and l ouf denote functions with momenta directed 
towards or away from the interface. The scattering matrix for holes is written 
in terms of the scattering matrix for electrons as S_ij(p\\) = Sji(—p\\) T . For 
negative Matsubara frequencies the boundary conditions coincide with the 
ones for advanced functions, which here simply means interchanging S «-> 
and the side indices (i «-> j). 

For definiteness, we consider a simplified spin-active model in which the 
scattering matrix is independent of the sign of the parallel momentum pu. 
In this case e.g. the scattering matrix for transmission can be written aa^ 

5i2 = 5 2 i = e^ 2 [82i + 4i (A • <?)] exp [i(fi ■ a)ti/2] , (11) 

where we have an average transmission parameter S21 = (\t+\ + |t_|)/2, 
a spin-filtering transmission parameter s 2 i = (l^+l — |i | )/2, plus a spin- 
mixing angle The amplitudes t± denote the transmission amplitudes of 
the interface for the two spin-directions in the basis where jl ■ a is diagonal. 
The phase 99 does not play a role in the following. 

For small transmission and small spin mixing the energy gap A has a 
step function form, in which case fi i n retains its bulk form at the interface. 
Similarly, the incoming function on the F-side retains its bulk form, i.e. it 
vanishes /2,£ n = 0. The resulting outgoing amplitude at the interface on the 
ferromagnetic side can be written as /2,ou* = {A s + sgn(e n )ylt fi-a)ia y , where 

A s = Aq cos A t = iA sin >&, (12) 

and the prefactor is Aq = \t + \ 7rA/|e n |. These amplitudes are the initial 
conditions for the singlet and triplet components in the ferromagnet at the 
interface. The spatial dependence of all amplitudes can be expressed in 
terms of the effective coordinate xq = x/\ cos6 p \. We obtain the following 
explicit expressions in the ferromagnet 

(13) 



fs{tn,Xe) 



As cos I — 



■A ■ I Xe 

i A t cos a sin [ — 
Xj 
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e -k»*9 t (14) 

ft±(en,x e ) = c s sgn(e n ) A t e- knXe fl ± , (15) 

where k n = 2\e n \/vf. Here we introduced the components of fi parallel to 
the z-axis, cosaz, and perpendicular, fix, see Fig. [T] (a). The coefficient 
c s = [1 + sgn(p/ ■ x) sgn(e n )] /2 selects the correct sign of the momentum 
relative to the x-axis. 

We conclude at this stage: (i) The singlet component is purely real while 
the triplet components are purely imaginary. It follows that in the complex 
plane f s (—z*) = f s (z)*, ft(—z*) = —ft(z)* for each p. E.g., the retarded 
functions have the symmetries f^(—e) = f^(e)* and ft(—e) = —ft( € )*- ( n ) 
Spin-mixing (finite $) is crucial for triplet components to be induced at the 
interface, (iii) The components ft± are induced only when ft is misaligned 
with respect to the exchange field (here Jz). 

To investigate the symmetry properties we expand the correlation func- 
tions in partial waves. Since we have rotational invariance around the x-axis, 
and the correlation functions only depend on cos P through the projection 
of the Fermi velocity on the x-axis, we can expand in Legendre polynomials 
Pi (cos P ) and get 

oo 

f s (cos p , e n ,x) = ^2 f s (l,e n ,x)Pi(cosO p ), (16) 

1=0 

21 + 1 f 1 

f s (l,e n ,x) = — - — J ^d(cosOp)f s (cosO p ,e n ,x)Pi(cosOp), (17) 

and similarly for the triplet components. To proceed we have to model the 
dependencies of the tunneling probabilities and the spin-mixing angle on 
the trajectory angle and then evaluate the various harmonics. The simplest 
model is a tunnel cone model with constant transmission probabilities and 
spin-mixing angle within a range of trajectory angles, i.e. = \f 7 T+, 

\t-(0 p )\ = ^/TZ, and $(0 p ) constant simply called i?. In the following we 
assume a very wide tunnel cone (— > vr/2), but it is straightforward to obtain 
the expressions for a more narrow cone. After integrations over cos P we 
obtain the amplitudes 

f s (l) = A [sgn(e n )] l [cos$ReQi(z) - cos a sin tHmQ^)] , (18) 
ftzil) = iA [sgn(e n )] m [cos a sin i? ReQi(z) + cosi? ImQi(z)] , (19) 
ft±(l) = M [sgn(e n )] m [sin^Q / (^)]/2±, (20) 

where z = k+x = 2(|e n | + iJ)x/vf. The function Qi(z) for the first few 
I is given by Qo(z) = zT(— l,z)/2 (s-wave), Qi(z) = 3z 2 F(— 2, z)/2 (p- 
wave), and Q2(z) = (5/2) [3z 3 r(— 3, z)/2 — zT{— 1, z)/2] (d-wave), where 



ftz(tn,X e ) 



c s sgn(e n ) 



At cos a cos 



iAc sin 
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Table 2 

The pairing correlations in ballistic superconductor-ferromagnet systems with spin- 
active interface scattering can be identified with the four types listed in Table [TJ 



Eq. ((HI), singlet: 


even I: even parity, even in frequency 


type A 


odd I: odd parity, odd in frequency 


type B 


Eqs. (P])-P0j). triplets: 


odd I: odd parity, even frequency 


type C 


even I: even parity, odd frequency 


type D 



r(n, z) is the upper incomplete gamma-function. Note the different spatial 
dependences of the singlet and m = triplet on the one hand, and the 
perpendicular triplets on the other hand, as they enter in the argument 
of Q[. Also note that Qi(k n x) is purely real since k n x is real. The Pauli 
principle and the symmetry requirements it imposes is affirmed in Eqs. f)18[> - 
(|20p . We can also conclude that we have all four types of correlations, see 
Table [21 including equal-spin pairing correlations when p, is not parallel to 
J. 

In the region £ j <C x -C £o of the ferromagnet the various components 
show different decaying behaviors depending on whether the pairing corre- 
lations involve two spin bands (f s and ftz) or only one spin band (ft±)- The 
asymptotic form of the gamma function for \z\ 3> 1 is T(n,z) ~ z 1l ~ 1 e~ z , 
and we can obtain simplified expressions for x 3> £j. We focus on the lowest 
Matsubara frequency, for which 

QlWUtj « ^|e-^e-^, (21) 

where £o = Vf/27rT c . The singlet and triplet amplitudes with zero spin- 
projection have the forms 



fs(l) = fo 

fUD = ifo 



sin jx/tr) . cos (x/gj) " 
cos 17 t- h cos a sin if — 

cos(x/0) . Q sin(x/0) 
- cos 17 t- cos a sin if — 



~ x ^°, (22) 
~ x/ Z°, (23) 



where the prefactor is / = \/T + T- (A/T) (21 + 1) /2. Note that there is a vr/2 
phase shift between the oscillations of the two components f s and ft z - Also 
note that the trajectory resolved functions decay exponentially on the large 
coherence length scale £o 5 while the various harmonics decay as (x/£;j) -1 , 
before the exponential decay on the £o scale sets in at large distances; the 
difference comes from the average over trajectories when projecting out the 
various harmonics. On the other hand, for the perpendicular triplets the 
asymptotic behavior is not reached until i>^i an d is of the form 

ft±(l)^ifosh^^e- x ^fl ± . (24) 
x 
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Fig. 2. Pair correlation functions in the ferromagnet of a ballistic superconductor- 
ferromagnet junction. The first partial wave (I = 0) for the lowest Matsubara 
frequency, sq — irT, is shown, (a) For the singlet, f s , and triplet with zero spin- 
projection on the exchange field, f tz , we observe a l/x-decay and oscillations, both 
on the scale £j = Vf/2J. On the same scale the perpendicular triplet f t ± only varies 
slowly, (b) At large distances all components decay on the coherence length scale 
£o- However, the magnitudes of f s and ft z are considerably reduced compared with 
ft± before this region is reached. We have chosen £ = 100£j, a — 7r/4, $ = 7r/4, 
and we have normalized the functions by (A/T)y/T+Fl. 

We plot in Fig. [2] the I = component of the Green's functions in the 
ferromagnet, as given by Eqs. f)18[) - (|20|) for the lowest Matsubara frequency. 
The higher order partial waves look very similar and have similar amplitudes. 
As can be seen in the figure and from Eqs. (|22p -(|23 p . in the region £j <C 
x <C Co the correlation functions f s and ft z decay like 1/x and are rapidly 
reduced by a factor Cj/Co compared to f t ±. A similar decay for f t ± is absent 
in that region. For x S> Co au components continue to decay according to 
x~ 1 e~ x ^° , as can be inferred from Eqs. (|22p -([24p. 

In fact, the Pauli principle requires odd-frequency amplitudes to be 
present in any inhomogeneous superconducting state, not necessarily spin- 
polarized. For example, the case of a normal metal coupled to a supercon- 
ductor is easily obtained from the above calculations and we recover for the 
case of a usual tunnel barrier (no spin-active scattering) the usual proximity 
effect. 2 - There is only a singlet component of the Green's function and in 
the normal metal region it can be expressed in terms of partial waves as 

/.(Z)= T^[sgn( en )]'Q,(^V (25) 
\e n \ V v f J 

which for large distances from the interface x 3> Co decays as ~ 
(Co/x) exp(— x/Co)- Above we introduced T, the transparency of the in- 
terface within the tunnel cone (here for simplicity chosen very wide — > tt/2 
as above). We see that higher partial waves I > 1 are induced, as is always 
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the case in inhomogeneous systems. In accordance with the Pauli principle, 
all odd components (I = 1,3,..) are odd in frequency. Thus, both singlet 
entries (types A and B) in Tabled] are always present in the usual proximity 
effect, and in fact in any inhomogeneous singlet superconducting state. Anal- 
ogously, the two triplet entries (types C and D) in Table Q] are characteristic 
for any inhomogeneous triplet superconducting state. 

3.2. Half-metallic Ferromagnet 

In this section we analyze the various symmetry components, intro- 
duced in the previous sections, for the case of a superconductor cou- 
pled to a half- metallic ferromagnet. This case was previously studied on 
theoretical grounds, and a non-zero Josephson-effect was predicted for a 
superconductor/half-metal/superconductor junction.^ The effect has re- 
cently been confirmed experimentally.^ In our previous article we focused 
on the p-wave triplet amplitudes, but it should be noted that there is also an 
odd-frequency s-wave triplet amplitude present in the half-metal. In fact, 
within our model all symmetry classes discussed in the previous chapter and 
shown in Table [1] are present in the heterostructure. In the following we 
discuss the results for the clean case. We have performed calculations with 
impurity scattering ranging from the clean limit to the diffusive limit and 
have confirmed that both the even-frequency p-wave triplet and the odd- 
frequency s-wave triplet components are always present and important for 
the Josephson current through the half metalP^ 

The full scattering matrix for our model is a 3x3 matrix, where the three 
scattering channels are the two spin channels in the superconductor and one 
metallic spin channel in the half metal. The second spin channel in the half 
metal is insulating and is not participating actively in transport phenomena. 
As discussed in the introduction, quasiclassical theory is therefore applicable. 
To be specific, we parameterize the scattering matrix connecting incoming to 
outgoing waves in the superconductor (S) and half-metallic (F) sides asj^ 



j,S 









( r n e£* r n e^TT-<?H) 



d ne ^TT+f) \ 



-TT / 



V dne-^TT-f) d iT e-^T+f) 

(26) 

where the extra phases of the transmission amplitudes ensure the unitarity 
of the scattering matrix, and the magnitudes are given by = 1 — £**/2W, 
r u = 1 - §/2W, r n = r; T = -t^/2W, r n = 1 - (t? T + t 2 ^)/2W, 
dff = ttf/W, and <% = t n /W, with W = 1 + (t^ + i| T )/4 The five real 
parameters of the scattering matrix above are t^\, t^, i?, i?^ and 

In the following we derive analytic expressions for a well defined limiting 
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case: the case of weak transmissions in all channels and small spin mixing 
angle. In the above scattering matrix this means that t+t and £it are 
small. Thus, in this limit the scattering matrix assumes the form, 





/ 


ea 





t ne i(*TT+S) ^ 




S = 









t lt e<(*iT-S) 


(27) 




V 


t ne -<(*TT-f) 




- 1 y 





Before presenting the results, we comment on the parameters in the 
scattering matrix. Within our theory, the scattering matrix is a phenomeno- 
logical input that characterizes the scattering of quasiparticles near the Fermi 
surface on either side of the interface of the system. Microscopic details of the 
interfaces are irrelevant for the low-energy physics near the Fermi surfaces. 
Thus, all the parameters are assumed to be independent of energy, and only 
depend on the positions of the scattering momenta on the Fermi surface. For 
strong ferromagnets, J is not a small quantity, and consequently the scatter- 
ing matrix has spin- active scattering terms of the order of J/ep, with €f the 
Fermi energy. Whereas the spin-mixing angle is a robust property of any 
interface to a strong ferromagnet, the spin-flip parameters t\t, i?it require 
the breaking of spin-rotation symmetry around the quantization axis of the 
ferromagnet. This can be the result of various mechanisms, for example the 
presence of interface regions with misaligned spin (magnetic grain layers). 
The above mentioned parameters represent in this case averages over the 
grain configuration along the contact region of the interface for each given 

sample, t^e M = {t\Ie l it)., where i numbers the different grains. Thus, 
there are expected sample-to-sample variations in the spin-flip parameters 
unless the typical grain size exceeds the contact size. 

We now proceed with the derivation of analytic expressions for the vari- 
ous symmetry components of the pairing amplitudes. On the superconduct- 
ing side of an interface with a half metal, all four components listed in Table 
Q] are induced and can be calculated perturbatively in the small parameters 
iff and t^. In the following we denote U n = + |A| 2 . We obtain for 
type D: fg = m$e n A/2Q,l, for type C: fg = -in$A/2tt n . In an expansion 
in spherical harmonics the two components f£ and fg correspond to s-wave 
and p-wave triplets. The type B component only enters in second order in 
t?, and has the form ff = 7rt? 2 e n A/4f^. Finally, the renormalization of 
the type A singlet component, 5ff = —ir-d 2 reduces the leading 
component only in second order in "&. Consequently, to linear order in i? the 
self consistent singlet order parameter A is not affected and stays constant 

up to the interface. All components, <5//, ff, fg, and f£ decay into the 

(S) 

superconductor exponentially with a decay length given by £n = u//2fi n . 
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Next we consider the equal spin pairing amplitudes in a symmetric 
Josephson junction with a half metal of length L extending from —L/2 < 
x < L/2. On either side of the half metal are singlet superconductors hav- 
ing order parameters |A|e* Xl and |A|e* X2 . The equal spin pairing amplitudes 
in the half metal, /ff are proportional to the components ft z in the super- 
conductors. For the even-frequency triplet amplitude we obtain in the half 
metal, 



f$ = -i P n sgn(cos 9 P ) , u A (e ,xa + e 4 * 1 ) + ^ (e** 2 - e 4 * 1 

V smh %7 cosh^f 

and for the odd-frequency triplet amplitude, 



(28) 



/fi = * A. sgn( e „) -^£r(e™ + + - , (29) 



where k n = 2\e n \/vf, fx = |cos0 p |, and (5 n = -nt^t^^ A| \e n \/2(e\ + |A| 2 ). 
The superconducting phases are renormalized by the phase shifts during 
tunneling, and are given by Xi = Xi ~ $i,TT ~~ W1 ^h i = 1,2. Note 
that equal-spin amplitudes are only possible when both ^ and ■d ^ 0. 
Even and odd frequency (i.e. p-wave and s-wave) triplets at the interfaces 
(x = ±L/2) are comparable in size. 

To illustrate the different symmetry components, we show in Fig. [3] 
numerical results for a superconductor/half-metal/superconductor junction. 
For these results we do not expand the interface scattering matrix in 9, 
iff, tii as m ^q- d2ZD, but use the general scattering matrix (f26|) and solve 
the boundary conditions numerically.^ We allow for an angular variation 
with respect to the interface normal of both the spin mixing angle and the 
transmission amplitudes; for definiteness we assume a variation proportional 
to | cos 9 p \. We iterate the order parameter and the boundary conditions until 
self consistency is achieved. We discuss first results for a 7r-junction (xi = 0, 
X2 = tt) as this is the ground state of the Josephson junction.^ To quantify 
the spatial dependences of the different symmetry amplitudes, we define the 
functions 



F a p(x) = -T (Pi(cos6 p ) fap{e n ,x,9p)^fh{e n ), (30) 

" e„>Q 



where Pi(cos9 p ) is chosen to project out the s-wave (1=0) or p-wave (1=1) 
parts of the pairing amplitudes, and /^(e n ) = 7rA/f2 n is the propagator for a 
homogeneous singlet superconductor that is introduced here to ensure con- 
vergence at large e n (note that at the boundary points f a p(^n) decays weakly 
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Fig. 3. Spatial dependences of all components of the pairing amplitude for a 
symmetric superconductor/half- metal/superconductor Josephson 7r-junction with 
interfaces characterized by the scattering matrix Eq. (|27p . The interfaces are in- 
dicated by the dashed lines. The components were numerically iterated until self 
consistency with the singlet order parameter was achieved. The interface parame- 
ters for normal impact angle are ■& = 0.37T, — 1, = 0.7, i9ff = = 0, and 
the temperature is T — 0.2T C . 



~ l/e n ). As is shown in Fig. El on the superconducting side all possible sym- 
metry components listed in Table [T] are induced, of which the even-frequency 
p-wave triplet and the odd-frequency s-wave triplet penetrate into the half 
metal. The latter two components are of similar magnitude, however the 
p-wave triplet is larger in the center region of the junction than the s-wave 
triplet. We have also performed calculations including impurity scattering 
and have shown that both components are essential for the Josephson current 
in the entire range from the ballistic to the diffusive limit. Details will be 
presented in a forthcoming publication.^ Also seen in Fig. [3] is that there 
are different symmetries with respect to the spatial coordinate, that result 
from the symmetry of the Josephson junction. The components shown are 
for the ground state; with a finite Josephson current (x2 — Xi 0> 7r ) an 
amplitudes become complex and the spatial symmetries are lost, as can also 
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Fig. 4. The same as in Fig. G2 but for zero phase difference. 



be inferred from Eqs. (|28j) - ([29j) . Note that in Fig.O as consequence of higher 
order terms in the transmission parameters and the spin-mixing angle, there 
are equal-spin pairing correlations also in the superconducting regions, and 
there are noticeable inhomogeneous contributions to the singlet amplitudes 
as well. 

For comparison, we present in Fig. U] the pairing components for the 
zero-junction case (xi = X2 = 0). In this case, the spatial symmetries are 
opposite to that of the 7r-junction. We have verified that the ir junction has 
the lower free energy.^ 

We consider next the Josephson current through the junction. We as- 
sume identical interface parameters for both contacts, that may depend on 
the impact angle 9 P via a parameter [i = \ cos0 p \. We obtain for small $, 
iff, and the following expression for the Josephson current density, 



J 

Jc 



-J c sin(x2 
eNfVf\A\*irT 




(31) 
(32) 
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Fig. 5. (a) The normalized critical Josephson current density J c (T)/J c o as com- 
puted from Eq. (|32p as a function of temperature T . The temperature dependence 
of the gap is given by the well-known BCS formula, (b) The zero temperature value 
Jco = Jc{T = 0) is shown in units of eNfVfT c \t^\ 2 \t^\ 2, d 2 as function of junction 
length L. (c) The peak position of the curves in (a) as function of L. 



where (• • -) M = f d/i(. . .), and Nf denotes the density of states of the spin- 
up electrons at the Fermi level in the half metal. We point out several 
differences to the case of the Josephson effect through a normal metal: (i) 
The minus sign suggests that the 7r-junction is stable for all temperatures and 
parameters within the above approximation, (ii) The Josephson current is 
proportional to the square of the spin-mixing angle i9 and the squares of both 
the tunneling amplitudes and iff. Thus, the magnitude of the critical 
current is much more sensitive to the interface characteristics than in a usual 
Josephson junction; in particular, strong sample-to-sample fluctuations are 
expected as the magnitude is proportional to t^. (iii) The additional phases 
$j.T+$Tt m Xi,2 can lead to a shift of the equilibrium phase difference between 
the superconductors except when they are identical for both interfaces of 
the Josephson junction, (iv) The Josephson current density is a result of 
both the even-frequency p-wave triplet and the odd-frequency s-wave triplet 
amplitudes, and neither of them can be neglected. 

In Fig. [5] (a) we show J C (T) from Eq. (I32p normalized to its zero tem- 
perature value. As above we assumed a variation of 6, \t^\, and |rjt| with 
impact angle that is proportional to | cos 9 P \. As can be seen, there is a low- 
temperature anomaly in J C (T) that has been discussed previously.^ Here we 
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show that this anomaly is a robust feature that exists even in the limit of 
small transmission and small spin-mixing angle, and is independent of these 
material parameters. Thus, the dependence on the interface parameters can 
be divided out and J c (T)/J c (0) is a universal function, only dependent on 
L/£o- The appearance of the anomaly can be traced back to the different 
energy dependence of the pair amplitudes in the half metal compared with 
a normal metal, which results in the e^-factor in the numerator of Eq. (|32|) . 
The position of the peak maximum in the temperature dependence of the 
Josephson current is shown in Fig. [5] (c) as a function of junction length 
L. It scales at large L as T pea k/T c « 2.5£o/L. For small L it saturates at 
a finite temperature. We also show in Fig. [5] (b) the variation of the zero 
temperature Josephson current with junction length. It decays for large L, 
as is L" 1 exp(— L/£q). 



4. SUMMARY 

In conclusion, we have investigated the superconducting pairing correla- 
tions with unconventional symmetries at interfaces between superconductors 
and ferromagnets. We have demonstrated that in ballistic superconduc- 
tor /ferromagnet junctions spin-active interface scattering naturally leads to 
all possible symmetry components of pairing amplitudes compatible with the 
Pauli principle. We have also discussed the case of a junction with a half- 
metallic ferromagnet. In this case odd-frequency s-wave and even frequency 
p-wave components are of comparable magnitude and are essential for the 
Josephson current. This leads to a 7r-j unction, where the supercurrent is 
carried by spin-triplet Cooper pairs. We have shown that a low-temperature 
anomaly in the Josephson effect through a half-metal is of universal nature 
in the tunneling limit in ballistic systems. We have derived analytic expres- 
sions for all pairing components and for the Josephson current in the limit 
of small interface transmissions and small spin-mixing angle. 
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